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THE NONLINEAR INTERACTION OF INTENSE LASER 
PULSES IN PLASMAS 


I. Introduction 

The interaction of ultra-high power laser beams^ with plasmas is rich 

2 

in a variety of wave-particle phenomena. These phenomena become 
particularly interesting and involved when the laser power is high enough 
to cause the electron oscillation (quiver) velocity to become highly 

relativistic. Some of the interesting laser-plasma processes which are 

3-12 

discussed include: a) relativistic optical guiding of the laser beam, 

b) the excitation of coherent radiation at harmonics of the fundamental 

13—16 

laser frequency, c) the generation of large amplitude plasma waves 

(wakefields), d) frequency shifts induced in the laser pulse by plasma 

17 18 19 

waves, ’ e) frequency amplification using an ionization front and 

f) single particle acceleration in a laser pulse. 

In the following, a fully nonlinear ID model is developed which 

describes the self-consistent interaction of intense laser pulses with 

plasmas. By assuming a "quasi-static" cold fluid plasma response, a set of 

coupled nonlinear equations is derived for the vector potential of the 

radiation field and for the electrostatic potential of the plasma. The 

quasi-static approximation assumes that in a frame moving at the speed of 

light, the plasma fluid sees a nearly steady-state radiation field (after 

the transients have decayed away). The resulting nonlinear equations are 

used to examine various laser-plasma interaction phenomena. The important 

2 3 20 

issue of laser-plasma instabilities ’ ’ is not addressed in this paper. 
Instabilities will certainly limit the laser p - r~-.gation distance and, 
therefore, will have a profound effect on some he potential 
applications discussed in this paper. 

Manuscript approved September 20. 1989. 
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Relativistic optical guiding is a result of the relativistic 
quiver motion of the electrons by the laser field. Analysis^ - ^ has shown 
that as the laser power exceeds a critical threshold, diffraction can be 
overcome, resulting in optical guiding of the laser pulse. Previous 

analysis of relativistic guiding have included, for the most part, only 

3-12 

the transverse electron motion in the plasma response current. 

Relativistic guiding was believed to occur on a fast time scale (on order 
of the inverse laser frequency). However, the present analysis finds this 
not to be the case. In the following nonlinear analysis of relativistic 
optical guiding, the electron density response and longitudinal electron 
motion are included self-consistently. It is shown that for short laser 
pulses (pulse lengths less than a plasma wavelength), the combined effects 
of the plasma density response and the longitudinal motion significantly 
reduce the relativistic guiding effect. It is found that relativistic 
guiding occurs only for long pulses with slow rise times (greater than an 
inverse plasma frequency). 

As the quiver motion of the electrons in a linearly polarized laser 
field becomes highly relativistic, the plasma response current will develop 
harmonic components. The harmonic content of the response current density 
can lead to the excitation of coherent radiation at harmonics of the 
fundamental laser frequency. It is shown that harmonic generation is more 
effective for short laser pulses (pulse lengths less than a plasma 
wavelength) than for long pulses with slow rise times. 

An intense, short pulse laser (pulse lengths near the plasma 
wavelength) interacting with a plasma can generate large amplitude plasma 

wave wakefields. The excited Wakefields may be used to: i) accelerate a 

13—16 

trailing electron bunch (laser Wakefield acceleration), ii) optically 
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guide a trailing laser pulse and iii) enhance the coherent harmonic 
radiation generated by a trailing laser pulse. 

Three other laser-plasma interaction phenomena are briefly discussed. 

The first describes how frequency shifts are induced in a laser pulse by a 
17 18 

plasma wave, ’ the second discusses how an ionization front may be used 

19 

to amplify the frequency of a laser pulse and the third describes how a 
single particle may be accelerated by the passage of a laser pulse. 


II. Nonlinear Formulation 

The ID fields associated with the laser-plasma interaction can be 
described by the transverse vector and scalar potentials, A^(z,t) and 
#(z,t), respectively (see Fig. 1). In what follows, we use the Coulomb 
gauge, i.e., 7-A = 0 which implies A z = 0. The transverse polarization of 
the laser field is arbitrary. The normalized potentials satisfy the 
following equations, 



!?)* ■ k p % ^ ■ k p k a - /Y> 


(la) 


3 2 <k 


3z 


2-t 


(r' O’ 


(lb) 


2 2 

where a(z,t) = |e|A^/m o c , <fr(z,t) - |e|$/m o c , kp = Wp/c, (Op = 

2 1/2 

(4n|e| n o /m Q ) is the ambient plasma frequency, n(z,t) is the plasma 

density, n is the ambient plasma density, 6. = v./c = a/y is the 

2 2 - 1/2 

normalized transverse plasma fluid velocity and y = (1 - (3 - §^) = (1 

2 1/2 2 1/2 

+ a ) /(I - 3 Z ) is the relativistic mass factor. In obtaining the 

right-hand side of (la) we used the fact that the transverse canonical 
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momentum is invariant and prior to the laser pulse interaction the plasma 
is assumed stationary. 

The fluid quantities n, 3 z and y are assumed to satisfy the cold 
relativistic fluid equations which can be written in the form, 

fa . c (nfy . 0, (2a) 


d t - - <1/2)1 ' -2 ( c h * It ) 3 - 2 + T <> - & H- < 2b > 


at ■ c6 z ff * It s a 2/3t ’ (2c> 

2 

where d/dt = 3/3t + cfi^d/dz. The term containing a in Eqs. (2b) contains 
the ponderomotive force. Thermal effects may be neglected provided the 
electron quiver velocity is much greater than the electron thermal 
velocity. Also, the ions are assumed to be stationary. 

It proves convenient to perform an algebraic transformation from the 
laboratory frame independent space and time variables (z,t) to the 
independent variables (?.,t), where £ = z-c3 t t and T = t. Here P t = 
v t /c = 1 is the normalized transformation velocity. To transform Eqs. (1) 
and (2) from z,t to ?,,T variables we note that 3/3z = 3/3J, and 3/3t = 3/3x 
- c0 t 3/3f,. Using these transformations, Eqs. (1) and (2) become, 


( -2 3 2 2P t r -2 d Z ) n 

h 73 * — H3T - c Z2>- r T a - 


8t 2 » U o 1 
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(3d) 
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2 2 2 1 n 

where y = l/(l-0 t ). Equations (3a-d), together with y = (1 + a ) / 

2 1/2 

(1 - 3 z ) , form a complete set of fully nonlinear, relativistic, cold 

fluid equations which describe the ID laser-plasma interaction. The ID 
model is valid as long as the radiation spot size is large compared to the 
plasma wavelength, i.e., r g >> = 2n/kp. 

In the following analysis, we choose v = c, since v = c lies between 
the group velocity and phase velocity of the laser pulse. Furthermore, we 
take the laser pulse to propagate in the positive z direction and to cross 
the z = 0 plane at t = 0. The £ = 0 plane defines the front of the pulse 
which extends into the £ < 0 region. In the £,x frame we are, therefore, 
interested only in the region where £ < 0 since for ( > 0 we have a = 0, 
n = n 0 * 3 Z = 0 and y = 1. 


a) Quasi-Static Approximation 

Equations (3) can be greatly simplified by noting that in the speed of 

light frame (0=1) we expect that under certain conditions a quasi-static 

state will exist in the macroscopic plasma quantities, n, 0 z and y. That 

is, if the laser pulse is sufficiently short, the fields a and 4> which 

drive the plasma are expected to change little during a transit time of the 

plasma through the laser pulse. From (3a) we find that the envelope of a 

changes on a characteristic time T ~ 2y|n /n|(co/« )/« , where w is the 
° e 1 o 1 ' p' p’ 

laser frequency. Since we will henceforth assume that « >> <*>p, the 
radiation envelope changes on a time scale which is long compared to a 
plasma period. If the laser pulse duration, t^, is small compared to t , 
i.e. << t , then the quasi-static approximation is valid. In addition, 
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the validity of the ID model requires that the laser beam diffraction time 

2 

(transverse spreading time), = nr s /(Xc), be long compared to T g . This 

condition is satisfied as long as r >> X . More formally, the quasi- 

s p 

static approximation (neglecting the right-hand side of (3c,d)) implies 
o 

s h J■>«' | « < 4a > 


\ I? I I « i. 


(4b) 


In this case, Eqs. (3c,d) can be integrated to give 


n(l-|3 z ) = n Q , (5a) 

r(l-P z ) - ♦ = 1. (5b) 

Using (5a,b) together with the expression for y, the coupled field 
equations become 


(1 i_ c -2 L.) _ k 2 JL 

Ic as, a-cj 3 t " p i++’ 

3-4 + k 2 ( 1 + <|>) _2 ( l + <J>/2) <f> = k 2 (l+<f>) 2 a 2 /2, 

ac p p 

where the plasma quantities in terms of the fields are 
n/n Q = 1 + \ ((l + a 2 )/(l + *) 2 - l), 

Y = (l+a 2 + (1+*) 2 )/(2(1+*)), 

3 z = (l+a 2 - (l + « 2 )/(ua 2 + (1+*) 2 ). 


(6a) 

(6b) 

(7a) 

(7b) 

(7c) 
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The above expressions for the fields and plasma quantities are fully 
nonlinear and hold for arbitrary polarization of the laser pulse. 
Consistent with the quasi-static assumption, a number of important points 
can be made concerning intense laser pulse propagation in plasmas. 


Ill. Laser-Plasma Interaction Phenomena 
a) Relativistic Optical Guiding 

The nonlinear index of refraction of the laser beam within the plasma 
determines, among other things, the optical guiding properties of the 
plasma. For the purpose of the present discussion, we assume that the 
laser field, a, is given by 

a = a L a,T)e ik */2 + c.c., (8) 


where a L represents the complex amplitude and k is the wavenumber. The 
characteristic spatial variation in the laser envelope, |a L |, is assumed to 
be of order L and is long compared to the laser wavelength, X, i.e., 
3|a L |/3£ = |a L j/L « k|aj = 2rt|a L |/X. 

Using the representation in Eq. (8) we find from Eq. (6a) that the 
refractive index, h = ck/w, is 

n = 1 - \ (k p /k) 2 /(1 + * s ), (9) 


where <|> is the slow part of the scalar potential and k << k. In 
s p 

obtaining Eq. (9), we replaced (l/(l+$)) g with l/(l + <J> g ), which is valid as 
long as kp << k (this implies that |^| << |$ s |» where <|>^ is the rapidly 
varying part of 4>). There are two cases that can be considered, depending 
on the envelope scale length, L, compared to the plasma wavelength, X . 
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i) Long Rise Time Pulse 

In this limit, L >> X = 2rtc/w , the first term on the left-hand side 
P P 

2 1/2 

of (6b) can be neglected and <p can be approximated by $ s = (l+|a^| /2) -1. 

The refractive index becomes 

h = 1 - \ (X/X p ) 2 (l + |a L | 2 /2)~ 1/2 . (10) 

Although the present analysis is ID, we expect that for a slowly 
varying transverse laser profile the index of refraction will depend on the 
transverse coordinates through the laser amplitude |a^|. Since the actual 
transversely-bounded laser beam amplitude falls off transversely, i.e., 

81 a^ | / 3r < 0, so will the refractive index, i.e., 8t|/3r < 0. The negative 
transverse gradient of the refractive index can lead to optical guiding. 

Since the refractive index is a function of the laser amplitude, the 
condition for optical guiding places a lower limit on ja^|. It is well 
known that if the refractive index is of the form given by Eq. (10), a 
critical laser power necessary for relativistic optical guiding exists and 
is given by 6 ~ 9 P crit = 17.4(X p /X) 2 [GW]. 

ii) Short Pulse 

Consider the short laser pulse limit, L < X p . When 1| << 1, Eq. (6b) 
can be solved for an arbitrary laser field a(C)> 

k ° 

<f> = -f J a 2 (£')sin(k p a'-S))d*;', (11) 

£, 

where the boundary conditions <(>(£=0) - 3<K£=0)/3£; = 0 have been used. 

If the pulse envelope is given by a^ = a^ Q sin(ix£,/L) for -L < l, < 0 and 
a^ = 0 otherwise, we find that the scalar potential within the laser pulse 
is given by 
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4> = (a Lo /8) { 1_(k p ~ Ar2/l2 '> jkpCos(2n£/L) - (4n 2 /L 2 )cos(k p £.)]|,(12) 

2 

where terms of Older (X/X p ) << 1 have been neglected. For L << X p , this 

gives <f> s = g(£)> where g(£) = £ 2 -2(L/2it) 2 [ l-cos(2n£/L) ]. Notice 

2 

that for L << X , $ is maximum at ?; = L where $ = (a. k L/4) . Also 

p s s ' Lo p 

notice that even for |a Lo l > 1> the assumption that <t> s << 1 is still valid 
as long as L << X . The index of refraction in the short pulse limit is, 
therefore, 

n = 1 - \ (X/X p ) 2 /[1 i- (Ji/2) 2 (a Lo /X p ) 2 g(*)], (13) 

where -L < £ < 0. 

In the short pulse limit, the fact that <f> s << 1 implies that the 

optical guiding effect is reduced significantly, by more than the factor 
o 2 

(ft /2)(L/X p ) << 1. The critical power, therefore, is increased by the 

inverse of this factor and, in addition, the degree of guiding varies along 
the pulse. Hence, it is unlikely that relativistic optical guiding can be 
effectively utilized in short, L < X p , laser pulses. 

Although it may appear that a long laser pulse may undergo guiding, 
assuming the various laser-plasma instabilities can be controlled, the 
front of the pulse will diffract. Initially, that portion of the head of a 
long rise time pulse in which the local power is less than will 

diffract. Once this portion has diffracted away, the pulse will exhibit 
"short pulse" diffractive behavior, i.e., the front region (~ X p ) will 
continue to diffract. The erosion of the front of the pulse due to 
diffraction will propagate back through the body of the pulse. The erosion 
velocity back through the body of the pulse (in the £ = z-ct frame) may be 
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2 

estimated by v £ = (X p /z R )c, where z R = nr g /X is the vacuum Rayleigh length, 
b) Harmonic Excitation 

The nonlinearities associated with the plasma waves can provide a 
source for the generation of coherent radiation at harmonics of the laser 
frequency. To examine this process we generalize Eq. (8) and represent the 
full radiation field by 

a = XI a^(£,x)e itIcJ 72 + c.c., (14) 

t 

where i = 1,2,3,... and a^ = a L (0 is the envelope of the dominant 
fundamental laser pulse, |a L | >> ja^|, for t > 2. 

It is clear from the right-hand side of (6a) that harmonic excitation 
is solely due to the excitation of the fast part of <j>. In particular, 
since the fundamental component of the radiation field dominates, the fast 
part of the scalar potential is 

4> f = - (k p /4k) 2 aL(l + <y“ 2 cos(2kO, (15) 

where we have taken the fundamental laser pulse to be a L (E,)cosk£. The 
source term in Eq. (6a) becomes 

S = k 2 a L (l + 4. s -Mj> f )“ 1 cos(kO = k 2 a L (l + <f> s )" 1 cos(kOE(«cos(2kO) m ,(16) 

H ^ m 

where is given by (15), a = ((k p /4k)a L ) 2 (l + 4> g )and m = 0,1,2,... 

As an illustration we consider the excitation of third harmonic 

radiation (3w). Since |4>^| << |<j> g |, the third harmonic component of the 

source is 

S 3 = (k 2 /2)(k p /4k) 2 aJ(l + 4» s )' 4 cos(3kJ;), (17) 
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where we have used the m = 1 term in Eq. (16). Substituting (17) into the 
right-hand side of Eq. (6a) and solving for the third harmonic field we 
find that 


l a 3 1 = j 4a LV’ < 18 > 

where x is the laser-plasma interaction time. The ratio of the third 
harmonic power to the fundamental laser power is 

2 

V P 1 ' [<V4X p ) 3 (l.+ s )- 4 a^ pT ] , (19) 

10 2 

where = 2.15 x 10 (r s a L /X) [Watts] (for a Gaussian transverse laser 
profile). 

Equation (16) shows that the generation of harmonics is a strong 
function of the plasma wakefield, described by 4> s , in the region of the 
fundamental laser pulse. For a single long pulse, large amplitude laser, 
|a L | >> 1, the slow part of the scalar potential is <f> s ~ |a L |//2. In this 

case, the harmonic content of the source term in Eq. (16) is exceedingly 

3 2 2 

small, and for the third harmonic, we have, P^/P-^ = (4(X/4A ) w^x/ l a L | ) • 

2 

Taking x to be a diffraction time, x^ = Jtr s /Xc and r g > X^, the third 

A 4 i i A 

harmonic power becomes P^/P-^ ~ (n/2) (X/X^) /(4|a^| ). In the case of a 

2 2 

short, L « Xp, large amplitude laser pulse, <j> s < |a L | (HL/2X ) , the third 

harmonic power is P^/P^ = ( |a^ |^(X/4Xp)^c0pX)^ = |a L |^(ix/4)^(X/Xp)^/4 for x 

2 

= x^ and r g > X^. For |a L | >> 1, a short pulse is more efficient than a 
long pulse for harmonic generation. 

Harmonic generation can be enhanced in regions of a plasma wave where 
4> s < 0. This can be achieved by positioning a short laser pulse at the 
appropriate location in a large amplitude plasma wave having a phase 
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velocity close to the speed of light. The large amplitude plasma wave can 

13—16 

be the wakefield generated by a iaser pulse or by an electron beam 

1 21 

pulse. 

c) Wakefield Generation 

The self-consistent evolution of the nonlinear plasma wave can be 
studied by numerically solving Eqs. (6a,b). Figures 2a,b show the plasma 
density variation 8n/n Q = n/n Q - 1 and the corresponding axial electric 
field E z for a laser pulse envelope given by a^ = sin( rt?,/L) for 
-L < ?, < 0. In this figure, L = X^ = 0.03 cm, X = 10 p, and a) a Lo = 0.5 

and b) a^ 0 = 2. The steepening of the electric field and the increase in 

22-24 

the period of the wakefield are apparent for the highly nonlinear 

situation shown in Fig. 2b (|a Lo | = 2) as compared to the slightly 
nonlinear case shown in Fig. 2a ( | a Lo I = 0.5). Figure 3 shows that the 
electrostatic potential <f> is predominantly slowly varying within the laser 
pulse even though 8n/n o has rapidly varying components. In addition, the 
figure shows that 4> can be negative behind the laser pulse. 

Qualitative aspects of nonlinear wakefield generation may be examined 
by considering a circularly polarized laser pulse with a square pulse 
profile, i.e., a^ = a Lo for -L < £ < 0 and = 0 otherwise [also, see 
Ref. 25]. For this case, Eq. (6b) may be solved analytically in terms of 
elliptic integrals. In particular, one can show that the optimal pulse 

length L = L (which corresponds to the maximum wakefield amplitude) is 

2 

given by L = 2y^ Q E(p)/kp •* 2y^ Q /kp for y^ Q >> 1, where E(p) is the 

2 2 2 2 

complete elliptic integral of the second kind, p = (y^ -l)/y^ o and y^ Q = 

2 2 2 
l+a^ Q . This gives a wakefield behind the pulse for which y^ Q > l + <j> s > 1 /y^q 

and the maximum axial electric field of this wakefield is given by E = 

6 max 

2 ~ 2 

(y^ o -1)/y 1o , where E = |e|E z /(m o c k ). The nonlinear wavelength of the 
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NL 2 2 

wakefield is given by X = 4y, E(p )/k -* 4y, /k for y, >> 1, where p = 

p lo o p lo p lo o 

4 4 

The large amplitude axial electric fields associated with the plasma 

waves can be utilized to accelerate an injected beam of electrons to high 
13—16 

energies (LWFA). In the region where Sn/n Q < 0, the transverse 

profile of the plasma wakefield can lead to a negative transverse gradient 
26 

of the refractive index, 3h/8r < 0. Equation (9) indicates that a 
properly phased trailing laser pulse (located at a maximum in <t> s ) may, 
therefore, be optically guided. From Eq. (16), it is seen that harmonic 
generation can be substantially enhanced in a properly phased (located at a 
minimum of <f> s ) short trailing laser pulse propagating in the wakefield 
generated by a leading laser pulse when -1 < <J> < 0. In addition, the sharp 

axial gradient in &n/n Q for a highly nonlinear plasma wakefield could 

17 18 

induce large frequency shifts in a short trailing laser pulse. ’ 


d) Laser Pulse Frequency Variations 

As the laser pulse propagates in the plasma, the excited plasma wave 

17 18 

can modify the laser frequency. ’ To analyze this effect we use the 
representation in Eq. (8) and set a L (E.,t) = |a^(£) |exp(i^(£, t)), where 
|a L (£)| denotes the laser pulse envelope, and \f»(£,t) is real and denotes 
the slowly varying phase shift. Substituting (8) into (6a) we find that 





T 

J (1 - ♦/(! + <f>))dt' . 


o 


( 20 ) 


The laser frequency is (^(S.t) = -9(k£ + vj»)/3t = ck + (c3/3£, - 3/3t)<j/ = 

2 

ck(l + (kp/k) /2) + Aw, where Aw is the frequency shift due to the plasma 
wave, 
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Aw = 


( 21 ) 


- T (k p /k) 2 (u* - c h I u* dT ')* 

o 

For completeness, we note that the wavenumber is given by k^(£, t) = 

3(k£ + ^)/3£, = k + Ak, where 

X 

* ■ I <V k)2 c ft I u| dT ’- < 22 > 

o 

and the refractive index, h = ck^/w^, is given by Eq. (9). The validity of 
(21) is based on the assumption that ^ is slowly varying; this breaks down 
when the frequency shift approaches the laser frequency, i.e., |Aw| ~ ck. 

For |<|>| << 1, <f> is given by Eq. (12) and the laser frequency shift for 
WpT >> 1 is given by 

o 

Aw = - ck(k p /k) 2 (w p x/4)k p | |a L U')| 2 cosk p (S' - f;)d£'. (23) 

Z. 

For a short pulse laser (L << X p ), assuming an interaction time x = x^, and 

spot size r > X , the frequency variation within the pulse is given by 
s - p 

3 2 2 

Aw/w = (n/2) (X(yX p )a^ o << 1. Frequency variations induced on a short 
pulse laser propagating in a plasma wave generated by a leading laser pulse 
can also be analyzed. ^ 

e) Frequency Amplification Using an Ionization Front 

A laser-induced ionization front may also be used to upshift the 
19 

frequency of a laser pulse. An intense pump laser with a short rise time 
(< 1.0 psec) may be used to ionize a gas as it propagates. This creates a 
moving ionization front in which the plasma density goes from zero to large 
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21 -3 

values (- 10 cm ) within a very short distance (< 0.01 cm). This sharp 

plasma density gradient propagates at the group velocity of the pump laser, 

2 2 1/2 

which may be approximated by v /c = (1 - w /« ) , where « is the 

go p 0 o 

frequency of the pump (ionizing) laser and is the plasma frequency of 
the ionized gas. One possible method for upshifting the frequency of a 
laser pulse is the following. A secondary short (< 0.1 psec) laser pulse 
propagates nearly parallel with the pump laser. By phasing the secondary 
pulse such that it "rides" the large density gradient of the ionization 

front, the frequency of the secondary pulse may be upshifted in much the 

17 18 

same manner in which a plasma wave may be used ’ to upshift the 

frequency of a laser pulse. For a linear gradient in the ionization front 

moving at velocity c, it can be shown that the frequency of the secondary 
laser pulse is given by 



where «(0) is the initial frequency of the secondary laser pulse, 

3(n/n Q )/3C = - 1/d, d is the spatial gradient length of the ionization 

front and z is the interaction length. As an example, consider a KrF laser 

21 -3 

(X(0) = 0.26 pm), an ionization density of n Q = 10 cm (X^ = 1.0 mm) 
with d = 0.01 cm and z = 1.0 m. This gives co(z=l m) = 26 co(0) or X = 0.01 
pm. Notice that it may be possible to tune the upshifted frequency by 
varying the interaction length. However, due to the difference between the 
laser pulse velocity and the ionization front velocity, the laser pulse 

will eventually slip out of the region of the ionization front, thus 

27 

limiting the interaction distance. 

19 

Alternatively, more dramatic "amplifications" in the frequency may 
be achieved by considering a secondary laser pulse propagating nearly 


15 






collinear to the pump laser, but in the opposite direction. In this case 
the ionization front of the pump laser appears as a relativistic plasma 
mirror moving towards the secondary pulse. Hence, part of the incident 
secondary pulse will reflect off this ionization front. The frequency of 

the reflected radiation, co g , will be relativistically doppler shifted, 

2 

w = 4(» /ft) ) ft)., where ft), is the frequency of the incident radiation, 
s o p 1 1 

This frequency amplification may be quite large, and one can envision c») g to 
be in the x-ray regime. The power of the reflected radiation, however, 
will most likely be small for two reasons: i) The power of the incident 
radiation needs to be sufficiently low so as not to appreciably ionize the 
gas, and ii) the ionization front may be a "poor" mirror, the reflectivity 
of which is dependent on the profile of the plasma density as well as on 
the ratio ft)./ft) . It should be noted that an electron beam may also be used 

l p 

in place of the pump laser to create an ionization front. 

f) Single Particle Acceleration 

As a final point we note that electrons, initially at rest, can be 
accelerated axially by experiencing only the transverse laser fields. The 
axial acceleration is due to the v x B force on the electrons. In the 
absence of a plasma, the scalar potential vanishes, <(> = 0. Equation (7b) 

shows that after experiencing the fields of a laser pulse, an electron, 

2 2 

initially at rest, will acquire a final energy given by = (a^/2)m Q c , 
where a^ is the final value of the laser's vector potential. This is 
reasonable since a^ is proportional to the area under the electric field. 
From Eq. (7c) and (7d) we find that the final ratio of the magnitude of the 
axial to transverse velocity is = (2+a^) /(2(l+a^) |a^ |) -* a^/2, 

where the second expression is valid for |a^ | >> 1. 
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IV. Discussion and Conclusions 


Based on a ID nonlinear quasi-static model, a number of different 
laser-plasma interaction phenomena have been discussed. Relativistic 
optical guiding is shown to depend strongly on the laser pulse duration. 

In the long pulse regime, optical guiding requires a minimum level of total 
laser power, i.e., the critical power. However, the leading portion of the 
pulse will experience diffraction. In the short pulse regime, relativistic 
guiding effects are greatly diminished by the density response and the 
longitudinal motion of the electrons. As the electron quiver motion 
becomes highly relativistic, the plasma response current develops 
harmonics. This process is analyzed and it is shown that coherent harmonic 
generation is more effective for short pulses than it is for long pulses 
with slow rise times. The generation of nonlinear plasma wakefields by 
intense, short pulses was examined, and wave steepening is observed for 
very intense laser pulses. Various applications of the plasma wakefields 
may be possible, including i) the acceleration of a trailing electron bunch 
(laser wakefield acceleration), ii) optical guiding of a trailing laser 
pulse and iii) enhancing the coherent harmonic radiation generated by a 
trailing laser pulse. In addition, the ID nonlinear model has been used to 
calculate the frequency shifts induced in a laser pulse by a plasma wave, 
as well as single particle acceleration by a laser pulse. The possibility 
of frequency amplification by using a relativistic ionization front may 
lead to a source of coherent x-rays. The viability of this process will be 
a subject of future research. 
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Fig. 2 (Con't) Density variation Sn/n^ = n/n Q _ 1 and axial electric field E z in 
GeV/m for a laser pulse located within the region -L < £, < 0, 
where L = = 0.03 cm, and a) a^ 0 = 0.5 and b) a^ 0 = 2.0. 
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Fig. 3. Electrostatic potential $ for a laser pulse located within the 
region -L < £ < 0, where L = = 0.03 cm and a Lo = 2.0. 
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